Adiabatically changing the phase-space density of a trapped Bose gas 
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We show that the degeneracy parameter of a trapped Bose gas can be changed adiabatically 
in a reversible way, both in the Boltzmann regime and in the degenerate Bose regime. We have 
performed measurements on spin-polarized atomic hydrogen in the Boltzmann regime demonstrating 
reversible changes of the degeneracy parameter (phase-space density) by more than a factor of two. 
This result is in perfect agreement with theory. By extending our theoretical analysis to the quantum 
degenerate regime we predict that, starting close enough to the Bose-Einstein phase transition, one 
can cross the transition by an adiabatic change of the trap shape. 



The observation of Bose-Einstein condensation (BEC) 
in magnetically trapped atomic vapors of rubidium jlj , 
sodium 0, and lithium Q has opened a new field of 
study at the intersection of atomic and condensed matter 
physics. Presently, condensates are produced routinely 
and detailed studies of condensate properties, such as 
collective excitations 0j and interaction of two conden- 
sates ||, are being made. The BEC phase transition 
itself is especially intriguing. Open questions include the 
kinetics of condensate formation and the effect of inter- 
atomic interactions and finite number of particles. Thus 
far, measurements have relied upon evaporative cooling 
to prepare the sample at the desired density below the 
critical temperature. Evaporative cooling, however, is 
inherently irreversible since it is based on the loss of 
hot particles from the trap 1^,0. To tackle the above 
questions it would be extremely valuable to vary the de- 
generacy of the trapped gas adiabatically in a reversible 
manner, with a fixed number of particles, especially since 
non-destructive detection methods have become available 

The possibility of increasing phase-space density and 
reaching BEC by changing the trapping potential was 
investigated by Ketterle and Pritchard For a colli- 
sionless gas they showed that it is impossible to influence 
phase-space density by manipulating the trapping poten- 
tial. An example is cooling by adiabatic expansion in a 
harmonic trap: one does not lose any atoms, but one 
does not get closer to BEC either. 

In this paper we show that this "no pain, no gain" 
principle is not true for a collisional gas. We experi- 
mentally demonstrate that the degeneracy parameter of 
a trapped gas can be changed adiabatically (without ex- 
change of heat) and reversibly, without sacrificing atoms, 
by changing the shape of the trap slowly compared to the 
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internal equilibration time. Our experiments are done in 
the Boltzmann regime where the degeneracy parameter 
nA 3 coincides with phase-space density. Here n denotes 
the density of the gas at the minimum of the potential 
and A = {2i:h 2 /mT) 1 / 2 is the thermal de Broglie wave- 
length at temperature T (with m the atomic mass and 
Boltzmann's constant /cb = !)• We observed a change of 
nA 3 by a factor of two, which agrees well with a quantita- 
tive prediction based on statistical thermodynamics ap- 
plied to a trapped gas. We point out that, given suitable 
starting conditions, this approach is also suited to cross 
the BEC phaseline. This follows from analytical expres- 
sions for heat capacity, entropy and condensate fraction 
of a Bose gas as a function of the trap shape. 

Our experiment is performed with atomic hydrogen in 
the cryogenic Ioffe trap described by Van Roijen et al jl0| . 
To determine quantities like temperature and density we 
measure the Lyman-a absorption spectrum of the gas 
and fit calculated spectra to the experimental one (see 
flTf for details). It takes a measuring time of 40 seconds 
to reach a 10% level of accuracy under present condi- 
tions. To minimize heating by photon recoil, we used our 
Lyman-a source at low intensity (typically 10 6 photons 
at the sample per pulse, 20 pulses per second). In order 
to assure sufficient signal to noise ratio, we replaced the 
photodiode used in previous experiments [Tl| ] by a pho- 
tomultiplicr. This method offers a non-destructive way 
to follow the evolution of a trapped gas in situ under 
changing conditions. 

The character of the Ioffe trap can be controlled via 
the magnetic field at the trap center, Bo. For /is -Bo 3> T 
(with [Ib the Bohr magneton), the trap is essentially har- 
monic. For hbBq <C T, the trap is close to linear in the 
radial direction and nearly harmonic axially. We selected 
two magnetic field configurations, A and B, shown in 
Fig. ^, that maximize the change in phase-space density 
within the constraints of our apparatus. 

After loading the trap from a cryogenic dissociator the 
sample is evaporatively cooled by ramping down the mag- 
netic barrier at one of the longitudinal ends of the trap 
(at z — 5 cm, see Fig. |l|). Evaporation is then stopped 
by raising this barrier. Hereafter we let the sample decay 
and equilibrate for about 1000 s to prepare a thermal 
sample that is better than 95% doubly spin polarized 
with typical densities around 10 11 atoms/cm 3 pc| . This 
density was selected to assure negligible sample loss dur- 
ing measurement of the spectrum. We gradually change 
from trap A to trap B or vice versa in 60 s which is 
much slower than the average elastic collision time of 5 s 
at this density. After each change of trap a spectrum 
is measured to determine the number of atoms and the 
degeneracy parameter. 

Changing the trap configuration from A to B, the 
temperature increased reversibly from Ta ~ 46mK to 
Tb ~ 60mK. Both Ta and Tb showed an upward drift 
of ~ lOmK to the final values mentioned above reached 
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after 5000 s. In Fig. || we plot the degeneracy param- 
eter nA 3 versus time for a sample cycled between trap 
A and trap B. It can be seen that the phase-space den- 
sity differs by a constant factor 2.05 ± 0.13 between trap 
A and trap B. The number of trapped atoms N versus 
time is plotted in the lower graph of Fig. |[ Although 
T and n differ considerably between the two traps, the 
measured value for N is seen to be conserved within 
experimental error. The slow decay of N was found 
to be described by N(t) = N(0)/[1 + N(0)Gt] as one 
would expect for a second-order decay process. The rate 
G = (1.34 ± 0.05) x 10~ 14 per atom per second is consis- 
tent with dipolar relaxation ]ic|jl^| . 

Our results can be understood within a dilute gas 
model for N atoms at temperature T trapped in a deep 
external potential, so that evaporation is negligible. Col- 
lisions keep the gas in internal thermal equilibrium. The 
number of particles in the gas is sufficient to enable a 
thermodynamic description. Since we change the trap 
potential slowly compared to the thermalization time, 
thermodynamic processes proceed reversibly. Since there 
is no exchange of heat or particles with the environment, 
the thermodynamic entropy S of the gas is constant. 
In our experiment, the influence of quantum statistics 
is small and interactions between atoms do not influ- 
ence thermodynamics as the mean-field interaction en- 
ergy is much smaller than T. The degeneracy parameter 
can be expressed in terms of the single-particle parti- 
tion function Z\ and the total number of particles N as 
nA*=N/Z 1 §. 

The internal energy E can be calculated from the par- 
tition function to give 

E=(^+J)NT, (1) 

where 7 = (T/V e )dV e /dT, with V e = N/n = Z X A 3 the 
effective volume. Here 7T is the average potential energy 
per particle. For many trapping potentials 7 is a constant 
independent of T and V e scales like T 7 . For example, 
for a box 7 = 0, for a harmonic trap 7 = 3/2, and for a 
spherical quadrupole trap 7 = 3. The canonical partition 
function Z± /N\ can be written as the exponential of the 



Helmholtz free energy E — TS (see, e.g., 13 1) to arrive at 



the following expression for the degeneracy parameter: 

nA 3 = exp f | + 7 - j^j . (2) 

Since in an adiabatic process S and N are constant, the 
phase-space density does not change unless 7 changes. 
It immediately follows that many of the most obvious 
ways of changing the trapping potential, like isotropic or 
anisotropic scaling of a harmonic trap, do not influence 
phase-space density. If we change 7 adiabatically (e.g., 
experimentally, by changing Bq in a Ioffe trap) , nA 3 will 
change as e 7 . In an extreme case, by slowly changing 
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the trap shape from square well to spherical quadrupole 
the phase-space density in the center of the trap can be 
increased by as much as a factor e 3 « 20. 

For the two traps used in our experiment the value 
of 7 is weakly temperature dependent. For the measured 
temperatures we find 7a ~ 1-79 for trap A and 7b ~ 2.53 
for trap B, the difference A7 = 7b — 7a after changing 
trap shape always being 0.74±0.01. This implies a change 
in phase-space density by a factor exp[A7] = 2.10 ± 0.02 
in perfect agreement with the measured value of 2.05 ± 
0.13. 

We emphasize that it should be possible to achieve 
considerable changes in the degeneracy parameter also 
around the BEC transition. For this purpose we extend 
our theoretical consideration to the case of a noninter- 
acting Bose gas in the degenerate regime. This is a good 
approximation also for a weakly interacting Bose gas as 
long as the gas parameter na 3 <C 1, where a is the scat- 
tering length. Below the critical temperature we have the 
additional condition that HqU <C T, with U = Awh 2 a/m 
the scattering strength and no is the condensate density. 
Expressions obtained for the entropy of an ideal Bose gas 
are still good approximations around and above the criti- 
cal temperature of a non-ideal Bose gas. Interactions will 
change the shape of the condensate, but its entropy will 
always be zero. Because the influence of the interactions 
on the entropy of the above-condensate particles is negli- 
gible in the binary collision regime the interactions only 
become important when the condensate fraction becomes 
so large that the mean field energy of the condensate 
changes the effective potential for the above-condensate 
particles. Gases used in current investigations, such as 
Rb, Li, Na and H are sufficiently close to this ideal gas 
limit that thermodynamics are essentially unaffected by 
interactions, unless the condensate fraction becomes ap- 
preciable. 

In the degenerate regime we continue to assume qua- 
siclassical motion of the atoms. For clarity we restrict 
ourselves to the case of a power-law potential, although 
our expressions can easily be generalized to include the 
Ioffe trap. A power-law trap is characterized by a density 
of states of the type p(e) = Ae 1 / 2 ^ (See @ and (l| for 
details). The scaling parameter A determines the size of 
the trap. E.g. for a harmonic trap A = ^(hui)~ 3 , where 
lu is the trap frequency. The parameter S governs the 
shape of the trap. It can be shown that in the Boltz- 
mann regime 7 = 6 independent of T for all power-law 
traps. 

The internal energy of a Bose gas in a power-law trap 
above and below the critical temperature T c M is given 
by 

E = NT(l + s) g -^f\ ;T>T C , (3a) 
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£ = W r|^?±£™m l+i ;T<T„ ,3b, 



2 7 9J +J (1) VI", 
where the Bose-Einstein integrals are expressed in poly- 



logarithms g a {x) = J2i=i x l °- The fugacity z = 
exp(^/T) with fi the chemical potential fT^ ]. Note that 
as we have T c 3/2+<5 oc AT, Eq. © is in fact independent 
of the number of atoms. Above T c the fugacity is given 
implicitly by 



N = AT^ s T(-+5)gs +s (z), (4) 



with T(x) the Euler gamma function. Below T c , z = 1 
and the number of atoms in the condensate, No, is given 
by N /N = 1 - (T/T c ) 3 / 2+s . In the high temperature 
(Boltzmann) limit the parameter nA 3 introduced earlier 
reduces to the fugacity z. The specific heat at constant 
particle number and constant trap potential can now be 
found by taking the derivative of E with respect to tem- 
perature. We obtain 

C = n(^+s) (f i+s (z)-f i+s (z)); T>T C , (5a) 

g \ / 'J 1 \ 3/2+<5 



C = Nl -+6\f i+6 (l)l-j ; T<T C , (5b) 

where we have introduced f K {z) = ng K (z)/ g K -i(z). For 
S = 0, f and 3 we find f 5/2 +s (1) = 1-284, 3.602 and 5.346, 
respectively. In Eq. (|5a| ) the second term gives rise to a 
discontinuity in the heat capacity at T c as already found 
by Bagnato et al. [^5), which appears only for S > 1/2 
because for S < 1/2 the function gi/2+5 diverges, and 
hence /3/2+a tends to zero, as T — > T c . This jump in the 
heat capacity was recently observed by the JILA group 
@. The entropy S obeys cZS/cff = C/T, and hence has 
a kink at T c for (5 > 1/2. We find that 



S = N(f i+B (z)-H*j)i T > T c, (6a) 

S = Nf i+s (l) [^y + ■ T<T C . (6b) 

For the homogeneous case (6 = 0) these expressions can 
be found in standard statistical mechanics textbooks (e.g. 
). The entropy as a function of T/T c is plotted in Fig. 



It is noteworthy that the ideal Bose gas in a power- 
law trap in three space dimensions is isomorphic to the 
uniform Bose gas in 25 + 3 dimensions. There is no heat 
capacity jump when the effective dimension is less than 
four. Remarkably, the processes which change the de- 
generacy parameter are just those which correspond to 
adiabatic changes of dimensionality! 

If one moves along an isentropic line in Fig. |^ by in- 
creasing S, one can take a dilute gas close to BEC through 
the transition. For example, by varying i5 from to 3 it 
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is possible to Bose condense a gas that had initially a 
temperature 12 times higher than the critical temper- 
ature for that trap and given N. Eqs. (JsJ) allow one 
to calculate the fraction of condensate particles Nq/N. 
Starting with a Bose gas in an ideal harmonic potential 
at T = T c , and reversibly changing the trapping potential 
to an ideal spherical-quadrupole, we would arrive below 
T c with a condensate fraction of 0.33. For the Ioffe trap, 
starting at T — T c at the limit of high we expect 
a condensate fraction of 0.25 in the trap at the limit of 
low Bq. Of course in these case T/T c may already be so 
low that the influence of interactions can no longer be 
neglected. 

Our method allows one to control the degeneracy of 
a Bose gas and gently pull it across T c and back. Cy- 
cling times are only limited by the elastic collision rate 
which for the alkali atoms has been demonstrated to be 
much faster than for hydrogen. Combined with non de- 
structive diagnostics this yields the unique possibility 
to study both condensate formation and destruction and 
to establish the presence or absence of asymmetries and 
hysteresis. 
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which is a subsidiary of the Nederlandse Organisatie voor 
Wetenschappelijk Onderzoek (NWO). M.W. acknowl- 
edges a TMR grant by the European Commission. The 
research of M.W.R. is supported by the Royal Nether- 
lands Academy of Arts and Sciences (KNAW). 



FIG. 1. Magnetic field profile of trap A and B radially cut 
through the field minimum (left) and cut along the principal 
axis of the trap (right). The horizontal lines indicate the 
measured thermal energy T. Trap A has a depth of 0.72 K 
and a field minimum Bo of 226 mT, Trap B has a depth of 
0.85 K and a field minimum Bo of 12 mT. 



FIG. 2. Experimentally determined phase-space densities 
of trapped atomic hydrogen (upper graph), and number of 
atoms -/V in the trap (lower graph) as a function of time, 
while the trapping potential was alternated between trap A 
and B. The curves are fits to a second order decay, with the 
constant ratio 2.05 ± 0.13 of phase-space density between the 
two traps. 



FIG. 3. Entropy of an ideal Bose gas in power-law traps 
as a function of T/T c (right) and, below T c , also as a function 
of the above condensate fraction (left). 
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